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Abstract 



The essence of the Thomas-Fermi model is the assumption that the local 
behavior of a many-body system can be approximated by that of a homoge- 
neous system. In this paper, we present the natural extension of the static 
Thomas- Fermi treatment of dilute Bose condensates, by describing the dy- 
namical behavior of the condensate in the same quasi-homogeneous approxi- 
mation. In particular, we calculate the dynamical structure factor 5'(q, lo) of 
a low-temperature condensate, confined in a harmonic oscillator trap. The 
result is a remarkably simple analytical function, which, with the proper inter- 
pretation, gives a powerful and insightful description of the scattering prop- 
erties of the BEC-system. 
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I. INTRODUCTION 



The achievement [|T|- of the long-elusive goal of atomic Bose-Einstein condensation 
has set the stage for exciting applications, such as tests of fundamental mean-field theories 
0] and the atom laser 0. The engineering problems associated with most applications 
are generally less severe if the condensate is highly populated and it might be necessary to 
create 'larger' condensates before a practical use of atomic-trap BEC can be demonstrated. 
Regarding this point, we note that recent experiments (see for example Ref. [0) report a 
significant increase in number of condensed atoms. Thus, the 'large condensate' or 'Thomas- 
Fermi' limit, aside from its convenient simplicity, is also of great importance from a practical 
point of view. 

The essence of the Thomas- Fermi model is the 'quasi-homogeneous picture', which is the 
assumption that the behavior of a condensate near a position r, can be approximated by 
that of a homogeneous BEC with chemical potential equal to the local effective chemical 
potential: 

^^{v)=^^T-V{v), (1) 

where /iy is the chemical potential of the trapped condensate and V"(r) the trapping po- 
tential. Although many authors have explored the physics and the limits of this 
description, most efforts have been limited to the static description of a dilute BEC. In 
fact, the few studies that have calculated dynamical properties H^l starting from a 



Thomas-Fermi description, have gone beyond the quasi-homogeneous picture in order to ob- 
tain a discrete spectrum. This allowed the description of interesting effects due to the finite 
size of the system, but caused the formalism to lose some of the simplicity of the straight- 
forward quasi-homogeneous description. Perhaps, the inability of the quasi-homogeneous 
description to calculate such finite-size effects is to be blamed for the lack in the literature 
of a simple quasi-homogeneous treatment of the dynamical BEC structure. 

In this paper, we present such treatment, and investigate the limits and restrictions 
inherent to the model. The simplicity of the results, combined with a proper understanding 



of the limits, gives a powerful formalism that can be used to interpret experimental data, as 
well as to check certain limits of more complicated computational schemes. Furthermore, the 
insight gained from studying the limits of validity, indicates how experimental data taken 
from a finite-size system can be interpreted in terms of the homogeneous system. 

The quantity that we choose to calculate is the dynamical structure factor, which rep- 
resents the information content of non-resonant scattering data about the dynamical many- 
body structure of the scattering system. More relevant to the atomic-trap BEC-systems, 
resonant light scattering gives a cross-section which in the off-resonant limit is also propor- 



tional to the dynamical structure factor [|19|, [|20[, ||2T|. More precisely, the single-scattering 



differential cross section d?a /(Elduo, where dVL is an infinitesimal solid angle and uj the energy 
transfer (/i = 1 in our units), is equal to 

where q is the momentum transfer, S'(q, uj) the dynamical structure factor of the many-body 
scattering system, and /(q) the scattering length that describes the scattering of an incident 
particle by an individual target particle. For non-resonant scattering, the scattering length 
can depend explicitly on the momentum transfer, whereas for off-resonant light scattering 
the scattering length is the large-detuning limit of the usual single-atom resonant scattering 
length (~ A7/A, where A is the inverse of the resonant wave number, 7 the width of the 
excited resonant state and A the detuning of the incident light). 

The structure factor is the Fourier-transform of the density-density correlation function, 

5(q, uj) = (27r)-^ J d^x d^x'dt' exp [-iq ■ (x' - x) - ujt'] (p(x', t')p(x, 0)) , (3) 

where p represents the density operator and ( ) denotes the thermally averaged expecta- 
tion value. A substitution of the spatial integration variables in to sum and difference 
variables, R = [x + x'] /2 and r = x' — x, gives 

5(q,cj) = j d^Ra{R; q,uj) 
where cr(R; q, u) = {271)-^ J dh dt' exp [-i(q ■ r - ut')] (p(R + r/2, t')p(R - r/2, 0)) (4) 



is a dynamical structure factor density. In the quasi-homogeneous approximation, we replace 
the correct structure density by o"^'(R^''(q, its value for a homogeneous system of chemical 
potential /i(R): 

S{q_,uj) ^ STF(q,^) = J d^Ra^^^\q,uj) , (5) 
which leads to an analytical expression, as we show below. 



II. HOMOGENEOUS BEG IN THE BOGOLIUBOV APPROXIMATION 

In treating a homogeneous system, the natural choice for a single-particle basis is the set 
of plane wave states, labeled by their wave vector/momentum k. Describing the interparticle 
interaction in the shape-independent approximation by means of a scattering length a, the 
second-quantized Hamiltonian operator reads 

k k,k',q 

where we have included the chemical potential ^, X = ina/m, V denotes the macroscopic 
volume of the system, and c, represent the annihilation and creation operators. 

For the purpose of describing the static properties of a BEC of N atoms with a coherent 
condensate of an average of A^o atoms in the k = 0-mode, we can replace Ck=o and cj^^Q by 
V^iVo. Provided the depletion is low, {N — Nq)/N <^ 1, we can neglect terms that contain 
less than two factors of In this approximation, only valid for dilute {na^ <^ 1) and low 

temperature Bose condensates (fcsT < Xtiq, where uq is the condensate density, uq = Nq/V, 
T the temperature and ks the Boltzmann constant), the hamiltonian is phonon-like: 

H (-/xno + ^nl) + ^{k'^/2m - fx + noA)4ck + 



The free energy of the system, F, also contains the entropy, which we denote by S, and F 



is equal to ^ F = {H) - TS. 



To lowest order, consistent with the assumption of low depletion, we can approximate 
the expectation value of H by the first term of (|^ 

{H) ^ V (-/.no + ^nl) , (8) 

and determine the chemical potential by minimizing F with respect to no, giving 

H ^ noA . (9) 

To determine expectation values involving non-zero wavenumbers, we return to (0) with (^. 
In doing so, we also find it useful to introduce Hermitian operators 



k 



V2 



which represent the fluctuations of the density and velocity field |23]. The canonical boson 
commutator relations such as cj^, Ck' = 5k,k', are equivalent to the requirement that $k 
and Ilk are canonically conjugate: [$k, I^k'] = ^5k,k'- With the fluctuation operators the 
Hamiltonian takes on the following form: 

H = ^(A;V2m)^($k$-k + H^kHk) + /x$k$-k • (11) 
k ^ 



In obtaining (|TlD , we have neglected a constant term, unimportant in describing the behavior 
of the system. 

At this point, we invoke the quasi-particle concept and introduce quasi-particle fluctua- 
tion operators H'^, defined as in ( p!0D with the particle creation and annihilation opera- 
tors replaced by quasi-particle creation and annihilation operators. We cast the Bogoliubov 
transformation between particle and quasi-particle operators in terms of the fluctuation op- 
erators. It can be seen that this transformation cannot mix 11 and $ operators because of 
time reversal symmetry. Indeed, the density fluctuation expectation value of a time-reversed 
state is equal to the expectation value of the state, whereas the expectation value of the 
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velocity-field fluctation changes sign. The only transformation that preserves canonicity is 
then a simple 'scaling' transformation: 

Hk = ftk' , (12) 

which is the Bogoliubov transformation in fluctuation operator notation. We choose a real 
so that all $ and H-operators are Hermitian. In addition, the isotropical nature of the 
system suggests that the scaling parameter should only depend on the magnitude of the 
momentum, Ok = a^. The value of the scaling parameter is determined by minimizing the 
free energy with respect to a^. In computing the expectation value of the hamiltonian op- 
erator, the temperature average of the fluctuation correlations depend on the quasi-particle 
occupation numbers v^, (n|^n'_,^) = ($'_k'^k) = (1 + 2t'fc)/2: 



dF d 



, (13) 



dak dak 

although the occupation numbers cancel out in optimizing and we find 

«fc = , (14) 

where Ek denotes the usual Bogoliubov dispersion relation, Ek = \J (A;^/2m -|- /i)^ — jj?. 
Unlike a^, expectation values of observable quantities do depend on the quasi-particle occu- 
pation numbers Vk = {exp[j3Ek] — 1)^^- In the last line of ([I^ ) we display the square of the 
scaling parameter, because is the quantity that appears in the expression for the dynami- 
cal structure factor, (p_qPq) no5q,0 + '^o('^-q'^'q) = '"'O'^q.O + ^0«g ('^'-q'^'q) • CoUSeqUeutly, 

we find for the dilute homogeneous EEC of low depletion, the following expression for the 
dynamical structure factor density: 

a^^-\c^,u;)^aln, {2n)-' J dt'exp(z^t')(<^>'-q(O<(0)) 

= no [(1 + u,)Sico - E,) + Sicu + E,)] , (15) 

Eg 



6 



which is the low-temperature generahzation of the well-known zero-temperature result p4 . 
Note that the dependence on the occupation numbers of ( p!5[ ) is reminiscent of stimulated 
and spontaneous photon emission. 



III. THOMAS-FERMI DYNAMICAL STRUCTURE FACTOR 

We calculate the structure factor S{q,uj) in the quasi-homogeneous approximation 
using the above derived result (|T5|) for the homogeneous BEC. In this manner we find 

= [ d'R no(R) [[1 + z/,(R)] 6 {to - E,(R)) - u,{R)5 {cu + E,iR))] , (16) 

where the position dependence of a;q(R) and i^q(R), defined as 



"'f^' = exp|3B,(R)|-l ■ 

stems from the R-dependence of the effective chemical potential /i(R). We remark that in 
the quasi-homogeneous description, the phonon-like (or collective mode-like) delta-peak in 
the spectrum ([16D , implies that the energy and momentum transfer, u and q, determine the 
spatial condensate region that is probed: the positions R for which -^^(R) = uj. 

We perform the spatial integration of Eq.(|l6|) for the simple example of a spherically 
symmetric harmonic oscillator trap with potential V{R), 

V{R) = Y^Lf , (18) 

where ujt is the trap frequency and L the extent of its single-particle ground-state, L = 
1 / y/mujT- The Thomas- Fermi expression for the condensate density in the low-depletion 
limit is easily obtained from the homogeneous result, Xuq ~ (P), which gives the following 
quasi- homogeneous expression for the condensate density ?t,o(R): 

„„m^mjfL^e<R^R„). (19) 
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The condensate radius Rq ( V^(-Ro) = I^t, or Rq = L'^2iJiT/ujT ), is determined from the 
condition that the integral of the density is equal to N, Rq = L(15aN/ LY^^ , where we 
have neglected the depletion |]2B|. Spherical symmetry reduces the integral (0) to a one- 
dimensional integral over the radial distance R. Finally, we substitute R by the effective 
chemical potential /i, 



R = RqJI- fi/fiT- 



(20) 



In carrying out this substitution, we replace no by /i/ A and the delta functions by 

6{uj ±E,)^6{i2- i2,{uj)) \dE,/di2\-' , (21) 



where fj,q{uj) is the effective chemical potential in the region where Eg or —Eg equals u, 



fig{u) 



q^/2m 2m 



(22) 



and dEg/dfi = q'^ /2m/ Eg. The resulting expression for the dynamical structure factor is 
simple: 

'l^(l)^Vl-^x(l + ,,^) ifc.>OandO</.,(o;)</.T 



STF{<i,Uj) = < 



1 (J]l_\ EiM. 

2L0T V o-L^ J 







( exp(/3|a;|)-l ) if ^ < and < figiuj) < /iT 

otherwise , 



which is the main result of this paper. In Fig. 1, we show the dynamical structure fac- 
tor as a function of energy transfer, for a fixed scattering angle or momentum transfer q 
corresponding to q^ /2m = 0.1/ir- The temperature of the condensate is different for each 
curve, ksT = (dot-dashed curve), 0.3/iT (dashed curve) and 0.5/iT (solid line). Notice 
that the ratio of the intensities distributed over positive and negative transfer energies, 
is sensitive to the temperature. Indeed, ( P^D satisfies the principle of detailed balance: 
S'(q, cu) = exp{Pu)S{—q, ~uj), which is a general result This temperature dependence 
implies that scattering experiments can directly measure the temperature of the atomic-trap 
condensates simply by comparing the scattered intensities at lo and —a; energy transfer, 

rn _ i^/kB 

- ln[5(q,a;)/5(-q -a;)] " 
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IV. DISCUSSION OF LIMITS OF VALIDITY 



Although the results of this paper are limited to dilute condensates of low depletion, 
the quasi- homogeneous description has a much broader range of validity. One necessary 
condition for the validity of the dynamical quasi-homogeneous description, is the validity of 
the static Thomas-Fermi model. In the low temperature region discussed in this paper, this 
amounts to the requirement that the size of the condensate exceeds the extent of the ground 
state, Rq » L, or equivalently, fiT » ^^^r ( 10' i0' Furthermore, approximating 

the structure density cr(R; q, cu) by the value of the homogeneous system is only sensible if 
the spatial variations of the BEG are imperceptible in the region probed by o"(R; q, cu), i.e. 
if q > l^^, where 1^ is the scale on which the condensate varies spatially. For example, for 
a harmonic oscillator trap, we could choose 1^ ~ -Ro/3 since yu(R) varies by approximately 
10 % from the middle of the trap to = -Ro/3. The resulting restriction, q > requires 
a scattering angle larger than the scattering angle for coherent scattering q < Rq^, so that 
the quasi-homogeneous model is only useful in describing incoherent scattering. Finally, we 
note that long-time fluctations, even on short distances, are affected by the finite size effects. 
Roughly speaking, if a localized perturbation creates an excitation that lives long enough to 
propagate to a region of different density and reflect back to the position of origin, then one 
should describe the fluctuations in terms of discrete harmonics or eigenstates of the finite 
system. The quasi-homogeneous description is valid either if the lifetimes of the excitations 
are sufficiently small so that the reflected excitation is damped out, or, if the scattering 
data probe the fluctuations only over time periods less than t^, the time needed for the 
reflected excitation to return. The restriction to short-time fluctuations can be achieved by 
using an energy resolution for the transfer energy that is less than or equal to Au ~ We 
estimate by assuming that the wave front of the excitation propagates at the local velocity 
of sound, c = \J ji/m, and we equate to the time necessary for the wave front to travel a 
distance /„, t^, ~ Iv/c. In the middle of the trap, c = cjt x 

/2o/V2, so that ~ and thus 
Aa; ~ t^^ ~ ujTy<{Ro/lv)- We conclude by stating that the above considerations indicate that 
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the Thomas-Fermi dynamical structure factor, as a function of the energy transfer, should 
be interpreted as a 'smooth' version of the measured structure factor and that we should 
compare intensities integrated over frequency intervals larger than or equal to Aa;. In reality, 
the estimate of Ao; depends on the region of the condensate that is probed, which in turn 
is determined by the value of the energy and momentum transfer. The middle of the trap. 



which is probed on the high-frequency side for |a;| (|a;| near y (q2/2m + /xt)^ ~ A*t ) requires 

the lowest value of Ao; (estimated above), whereas on the low frequency-side, ~ /2m, 
the edge of the condensate is probed, where the quasi-homogeneous description cannot be 
trusted and Aa; — > oo. 
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Figure Caption 



Figure 1: Plots of the Thomas- Fermi dynamical structure factor as a function of the en- 
ergy transfer cu, for fixed momentum transfer q, /2m = 0.1 fix- The three curves show 
the results at different temperatures, ksT — for the dot-dashed line, ksT — O.S/it for the 
dashed curve, and ksT — O.B/it for the plot shown in solid line. 
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